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I Introduction 



The statistical mechanics of entangled polymers, i.e., polymer chains under topo- 
logical constraints is a generally unsolved problem. The main difficulty is to specify 
distinct topological states of the polymer chain. Closed polymers, or polymer loops 
appear to be a much simpler system as they are either linked (with themselves or 
with one another) or unlinked. Linear chains, however, can always be disentangled. 
On a shorter time scale than the disentanglement time though, it seems justified to 
define topological states "on the average" , using the same formalism like for polymer 
loops. 

Mathematically, the problem of specifying topological states of polymer loops 
is equivalent to the classification problem for knots and links [^. Since the mid- 
eighties considerable progress has been made following Jones as various new knot 
polynomials have been discovered. (For a review on knots see [0]). 

For an analytical theory of the polymer entanglement problem, the algebraic 
form of these invariants is not suitable (see section V for new perspectives). They 
are generally expressed in one, two or three variables which appear in the defining 
relations (known as skein relations). There is no immediate relation of these variables 
to the polymer conformation, and consequently there is no reasonable way how to 
couple the knot invariant to a statistical weight for a given polymer conformation. 
Algebraic knot theory seems to be applicable for the theory of entangled directed 
polymers Q], and certainly does so in computer simulations (e.g. [§] [0] ) • 

The degrees of freedom appearing in the statistical weight of given polymer 
conformation are usually expressed in terms of segment positions r(s) which is a 
mapping [0, A^] — > from the contour variable s to d- dimensional space [§]. For 
(i = 3, one invariant showing an explicit dependance on these variables is the the so 
called Gauss invariant for two given closed loops Ca and Cp parametrized by ^^{s), 
r/3(s) 

*(C„, C) ^ ^ I * is'.Js) A . ^fff^ . (1) 

which is invariant with respect to continuous deformations of the loops, was ffist used 
by Edwards to discuss entangled polymer loops [0 . It is also called the Gaussian 
linking number (sometimes also the winding number). As Edwards already noted, 

2 



in (i = 3 the Gauss invariant does not uniquely specify a given link. There is, in 
fact, an infinite series of higher order link invariants that appear quite naturally 
in a perturbative expansion of a field-theoretical representation of the generalized 



Jones polynomial [|11[] [|12| . For special cases such as a random walk winding around 
an obstacle of infinite length, rigorous results can be obtained using the Gauss 
invariant, as it has been discussed in detail by Wiegel | |13[ |. 

As even the Gauss invariant is difficult to handle mathematically for a rigorous 
treatment of the three-dimensional entanglement problem, many mean field type 
arguments have been used ||^ for a rough characterization of the topological states. 
Tube models for entangled polymer melts are the most prominent approaches, which 
seem to lead to contradictory conclusions for open linear polymer chains and closed 
ring polymers. Therefore more detailed knowledge about the topological states 
of polymer systems is needed. In general, this seems to be a difficult task, and to 
overcome mathematical and conceptual difficulties simple and model type situations 
must be studied, to learn about more complex ones. Only recently more detailed 
situations are under consideration by using the path integral approach and the 
Gaussian linking number constraint ||14|| . 

To give an example that shows the complexity of the structure of the theory 
we mention the following development: The "easiest" topological arrangement of 
closed polymer rings is the non-concatenated melt of rings, since all winding numbers 
between different rings are zero. Scaling arguments for the typical size of a ring, 
R oc A^^, have been put forward [T^, giving an estimate of u = 2/5. This result 
seems to be in rough agreement with computer simulation [|16| (see also ||17||[p^). 
The more detailled analytical many chain theory in supported the scaling result, 
although the theory is more involved than simple scaling. 

An important step forward to formulate the problem by field theoretic methods 
was made by Brereton and Shah |T^. A test loop in the melt was considered 
that is entangled with many other chains. The resulting theory can be mapped to 
Euclidean electrodynamics in 3 = 2 + 1 dimensions coupled to a 0{n) 0^ theory for 
the conformation of a self-avoiding walk (in the limit n |^ 0)- It is the basis 
for work by Nechaev and Rostiashvili [p^ in two dimensions. 
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As a matter of fact, for d = 2 the Gauss invariant becomes rigorous for "simple", 
i.e. non-self-intersecting loops. It reads 

G-^{C) = dsr{s) ■ V (In \r{s) - r,|) A ry (2) 

zvr Jc 

Here the vector r(s) represents the segment positions of the polymer loop in the 
plane. In Eq.(|I]), the role of the polymers entangled with the loop is taken formally 
by obstacles at the positions rj. rj is a unit vector perpendicular to plane |^ . 
Gi{C) is also called the winding number of the loop. Eq.(0) may also be expressed 
in terms of a Cauchy integral in the complex plane and is related to the oriented 



area of the loop in the plane |23|. For so-called "complex" loops, i.e. loops with 
points of self-intersection, special care is needed as the invariant Eq.(^ might give 
zero although the loops are entangled (see Rostiashvili, Nechaev, and Vilgis p^). 

The theoretical basis of this paper is the field theory by Brereton et al. [|1^]. It 
was studied by Nechaev and Rostiashvili in order to discuss the behavior of 



a polymer loop in an array of randomly distributed parallel line obstacles p2| ||24 
whose spatial distribution is quenched. A first order transition for critical length 
Nc for the polymer loop was found when the quenched average over the winding 
number distribution was taken. It has been interpreted as a collapse transition (for 
> Nc) to an octopus conformation resembling a randomly branched polymer with 
no self- interaction with an end-to-end vector scaling of the size R oc A^^:, where A^ 



is the total length of the ring [^ [p7 |. 



In fact, in the work below it will be shown explicitly (within the approximation 
made in |^) that the conjectured scaling behavior is valid even for annealed disorder 
with respect to the spatial distribution of obstacles. This assumption is reasonable if 
the obstacles are supposed to represent other polymers entangled with a given loop. 



and was also made in the original theory by Brereton et al. |T9[. We will obtain the 
same free energy (in the mean-field approximation) as Nechaev et al.[^ - yet in a 
physically more transparent form - and consequently, the same transition behavior. 

We will exploit the fact that the interaction due to the topological constraint is 
a so-called "area law", i.e. it is proportional to the area enclosed by the loop. This 
is an exact result known in quantum field theory in the context of Wilson loops 
and the confinement problem for abelian gauge fields (see e.g. PSI). Therefore, the 



area of the loop is equal to first order to the effective potential found in [^] and is 
dependant on the conformational fields and topological quantities. 

The paper is organized as follows: In section II we clarify that the effective 
interaction is proportional to the area being the order parameter of the problem. In 
section III, we show on the level of polymer field theory how the area appears in the 
effective interaction and how it depends on the conformational fields (in the mean 
field approximation). The critical behavior of the order parameter as a function of 
topological parameters is discussed. A new upper bound for the range of stability 
for the mean field solution (in the replica-symmetric case) is found. In section IV, 
the conjectured scaling like for randomly branched polymers is shown, and in section 
V a brief outlook on the complete, i.e. d = 3 entanglement problem is given. 

II The area as the order parameter 

In the work of Nechaev and Rostiashvili [^] it is not evident at first sight that the 
first order phase transition caused by topological disorder corresponds to a classical 



collapse transition, similarly to the case of a polymer immersed in bad solvent ^1 
We expect indeed structural differences between a collapsed chain in bad solvent 
and a "collapsed ring" under topological constraints. 

In order to understand better the nature of the phase transition, it is desirable 
to introduce an order parameter that fits to the problem, which is in our case the 
area of the 2d projection of a simple loop. Following Cardy we use the covariant 
expression 

A = Jd^^J rfV(A^(x)A,(x')) J^(x)J,(x'). (3) 

where x = (x, y) is a vector in 2-dimensional Euclidean space. Variables are the 
polymer current densities or tangent vector densities Jfj,{'x.) = /J^ (isr^(s)5(x — r(s)) 
where N is the chemical length of the chain. The gauge fields are of the U(l) 
type and their correlator is gauge dependant. Using the gauge Ai = it reads as 
{A^{x.)A^{x')) = —■^6^o6uo\y — y'\S{x — x'), and it becomes obvious that the r.h.s of 
Eq.(^ is indeed an area. 

When using expression Eq. (RT) for the definition of the area some care is needed. 
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The above equation is valid only for non-self-intersecting loops which will be the 
scope of the present treatment. In the case of self- intersections, negative area con- 
tributions may cancel positive ones giving a total zero area (see above the discussion 
following Eq.dl)). 

We next recall the incorporation of constraint Eq.(^ in the partition sum for 
the loop For a fixed number of obstacles c enclosed by the loop the partition 
function is given by: 

Z{c) = ypr(s)5(r(iV) -r(0))5(^c-£dsr(s) ■ a) (4) 

exp (^~'j2 f (^^^"^i^) ~ ~2 '^'^'^(^('^) ~ '^(■^0)^ 

/ is the Kuhn segment length, and is the 2d excluded volume. From Eq.(^ we 
see that the gauge field A is given by I^ili V (In |r(s) — r^l) A ry so that the delta 
function fixes the winding number in terms of the 2D Gauss invariant. It satisfies 
V ■ A = and V A A = ry {(p{r) — (^q) where ip{r) = J2i <^(r — ^i) and ipo is the mean 
density of obstacles in the xy-plane. 



In contrast to |^ we suppose that first, the spatial distribution of obstacles 
is annealed whereas second, the distribution of winding numbers c, i.e. obstacles 
enclosed by the loop is quenched. The first assumption is reasonable if the obstacles 
are to represent other polymers (of infinite length) entangled with the loop, which 
is an important feature of a model that projects the original three-dimensional en- 
tanglement problem to two dimensions. 

The second assumption is clear from the fact that once a given winding number c 
is fixed for the loop, it should remain fixed in the process of averaging over both the 
conformations of the loop and the positions of obstacles. As a consequence, we take 
the annealed average over the spatial distribution of obstacles to be the gaussian 

P (^(r)) ~ exp (-^ / ^'x(^(r) - ^o)') ~ exp ("^ / d'^V A A)^) . (5) 

The quenched distribution of the number of obstacles c is assumed to be a gaussian 
with mean cq and dispersion Ac, 

P(c)~exp(-(£^) (6) 
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and is used to average the free energy. 

The winding number constraint in the partition sum is expressed by a Fourier 
transform introducing the variable g, a chemical potential conjugate to c: 

S{c- £ dsr{s) • a) = 1 ^e'^^'-'^fc '^-(^)-a (7) 

The distribution of winding numbers -P(c) may then be transformed into a distri- 
bution P{g) for the chemical potential p2[. For later purposes it is crucial to note 
that the g"^ averaged over P{g) is 

The partition function is now expressed in terms of g. After averaging over the the 
spatial distribution of obstacles, it reads: 

{Z{g))A = Ary'l?A5(V-A)|Pr(s)5(r(iV)-r(0)) 

exp I — [ (i^x(V A A)^ -ig (f dsr{s) ■ A 

\ 2(fo J Jc 



^ j dsi\s) -'^jdsj ds'5{v{s) - r(s')) j (9) 



TV is a normalization factor for the average over the gauge fields . 
Carrying out the integral over the gauge fields A, one obtains: 

{Z{g))A = Jvv{s)5{r{N)-r{0)) 

exp (^-^ J rf^x J rfV(A^(x)A,(x')) J,(x) J,(x') 

~ j dsr\s) -^f'^'^f ds'6{r{s) - r(s'))^ (10) 

The resulting term in the exponential is proportional to the area of the loop. In 
fact, the interaction reads as 

= (11) 

If we replace g^ by its mean value Eq.(H) we obtain: 



The approximation considered in |^ is limited to range of values for ipQ and 



which require the factor in front of A to be always positive (as to further restrictions 



on the set of values for these parameters where a mean field solution is valid, see 
the end of section III). Consequently, in order to minimize its energy the loop tends 
to collapse, decreasing its area. 



Ill The area as a collective variable and polymer 
field theory 

We now introduce the area explicitly as collective variable in the partition sum using 
1 = / dA6{A — A) where A is given by Eq.(|^). After transforming the delta function 
and some standard manipulations we then have: 

{Zig))j^ = ^^|dA|da| Vv{s)5{v{N) -riO)) 

exp (^taA-{ta + ^) J rf^x J d^^'{A^{^)A,{^'))J^{^)J,{^') 

-^j dsi\s) -jfds<f ds'6{r{s) - r(s'))^ (13) 

2 

Let us define ia = ia + ^^y- for a moment, and express a in terms of a, and make 
the change a — a in the integration variable. This procedure looks strange at first 
glance, because the new integration variable a becomes now complex. Below we 
show that this is not a serious problem for the purpose in this paper. After these 
manipulations we obtain: 

{Z{g))A = AT J dA J da J Vr{s)6{r{N) - r(0)) 

exp (^-^^A + i&A-ia J d^^ J rfV(A^(x)A^(x')) J^(x) J^(x') 

f dsr^s) -'^jdsj ds'5{r{s) - r(s')) j (14) 

In the next step, the term depending on the gauge field correlator is expressed in 
terms of a gaussian integration: 

{Z{g))A = ^ J ^^A5(V ■A)JdAjdaJ Vr{s)6{r{N) - r(0)) 

exp l^-i J d^^iV A A)2 -iej d^xA^{^)J^{^) - + zaA 

j dsr^s) -^f^^^f ds'5ir{s) - r{s'))\ (15) 



e is a shorthand notation for \/2ia and is the " couphng constant" of the analoguous 
Wilson loop problem, well known in quantum field theory p8[. We note again that 



e is complex, but it will turn out below that this is not a problem. The partition 



sum has now a structure similar to the original formulation considered in p2[ and 



is suitable for a field theoretic treatment. Note that the interaction ^^y-A has been 
completely separated from the conformational average and the average over the 
distribution of obstacles. The area is related to the conformation directly only via 
the coupling constant e or a respectively. In fact, the problem will be first examined 
for a given realization of these variables. In the last step, the parameter a will be 
eliminated to give back the dependance of the area on the conformation of the loop. 

To proceed further, we consider only consider the functional intregation over the 
positions r(s), i.e. the partition function 

Z(e; [A]) = AT y" Pr(s)5(r(iV) - r(0)) exp (^-^ ^ dsx'^is) 

-y/rfs fds'6{r{s)-r{s'))-te J rf2xA^(x) J^(x) j 

The partition function Eq.(^) describes the statistics of the loop for given external 
"magnetic" field A. It is formally the same partition function as Z{g) in Eq.(P), and 
can be treated in terms of the n vector 0^ theory in the limit n —* according 
to |T^[P^]. To obtain the field theory, the following standard steps have to be 



carried out. First, the two-dimensional excluded volume interaction is expressed in 
terms of a gaussian average over a pseudopotential Then, one has to consider 
the Green's function of the loop for a given realization of the gauge field A and the 
pseudopotential for the excluded volume. The Green's function is expressed in terms 
of a Gaussian field theory. The averages over the pseudopotential and the gauge 
field lead to consider an n-fold replicated field theory (see |]22| for technical details). 
When the average over the pseudopotential is carried out, one finally obtains 

Z„(e; [A]) = n (/ J exp J d^Kn[A, , (16) 



with 



(17) 



The fields 0j are replica fields for polymer loops. The model has been embedded into 
3-dimensional space, so the excluded volume term is the embedded 2-dimensional 
one with L being the mean size of the polymer in the z direction. 
Now, the average over the obstacle distribution is taken as follows: 

(Z„(e; [A]))a =M j VMiy ■ A)Z„(e; [A]) exp {-]^ J d^^{V A A)^) (18) 



The gauge fields are integrated out following |22]. To simplify the algebra, the 
Landau gauge is used. One obtains an effective action with the one-loop correction 
being given by: 

(19) 



^i-ioop = / ^ log (fc' + g E 



with Q = = iaP. As Q is complex, a complex logarithm in Eq.(|l^) has to 
be considered. It is easily shown that by restricting the analysis to one Rieman- 
nian sheet, the integral can be evaluated straight-forwardly under the assumption 
that the 0, are constant in space. This assumption is consistent with the result 
that the topologically restricted chain forms a dense object with very small density 
fluctuations. This corresponds indeed to the assumption of (pi ~ constant. 

In the replica-symmetric case is studied because only in this case the effective 
potential approximation can be used. (For the details of solving integral Eq.(|19[) 
and the renormalization procedure we refer the reader to [p^2|| ) . 



In this case, we now approximate J2i4'i4'i = ncfxp*. After renormalization ac- 
cording to p2[ one obtains: 

{Zn{e = e(«); [A]))a = exp [- j c/^^R^e//) (20) 
with an effective Lagrangian 

(21) 

M is an arbitrary subtraction point appearing due to the renormalization procedure. 
Because the fields (j) are now constant in space, the integration in Eg. (|21|) gives simply 
a constant volume factor V . From Eq. (pO|) one obtains the contribution to the free 
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energy as a function of a which is conjugate to the area. Devided by the system 
volume, it is given using the standard formula: 

/(") = ^ = ^|^(^n(e = e(«);[A])UUo 

(22) 

The next step is to transform back from a to the area A. This is done by a Legendre 
transfrom (or by a Fourier transform of the partition function). Finally one has to 

2 

add the area term ^^y-A to the free energy which yields the partition function 
averaged over the distribution of obstacles: 

(^(^))a = / dAb [a + ln(M!) _ ^\\^\^ e-^/(^'^) (23) 

with the free energy density 

= + {m- - La^M^m + (24) 

The set of Eq.s (^) and ( p4|) is the fundamental result of this paper. The free 
energy density f{A,g) is indeed the area law plus the renormalized action for the 



self-avoiding walk loop. Integrating over the area gives back the result of for 
the free energy density in terms of conformational fields and g only. Averaging the 
free energy density in Eg . (|2^) over the distribution of winding numbers using Eq.(|]) 
one finally obtains 

f{A) = IfiAg)], = ^ (l - ^ + - LaVMl^ + ^\<j>\^ (25) 

As a consequence of Eq.(^), the essential result we have obtained here is the de- 
pendance of the area on the fields in the mean field approximation expressed in 
the delta function of Eq. (|2^) . 



Introducing the segment density p = |0p and choosing = we then obtain: 

A = yi^(p-ipln(pL^)). (27) 
11 



Let us now investigate the behavior of the area when the collapse transition occurs. 



It has been studied in detail in ||2^, so we just mention the results. In fact, the 
collapse transition takes place at the critical length Nc which is given by (see ^2 
but with the factor replaced by ^): 



^ = ^^(l-^] In 



N, 8n A, 



A. 



4:71 A, 



A. 



As 1/Nc is always positive, the condition 



2 72 



-) ^^11-^1 >1 



47r A, 



A, 



(28) 



(29) 



follows for the stability of the mean field solution. Ai N = N(. the segment density 
is: 



Pc 



La'^ ATI A, 



A. 



Thus the critical area is: 

72 \ 2 



V 



2La? \ 27r 



^0 

A 



A, 



a J 47r A, 



A, 



(30) 



(31) 



Ac essentially depends on the topological parameters cq, the mean winding number, 
Ac, the dispersion of the winding number distribution, and y^o the mean density of 
obstacles in the plane. As A^. should remain non-negative we obtain a new upper 
bound in addition to the inequality Eq.(E^): 



Atx Ac 



Ac 



(32) 



where e= 2.714... is Euler's constant (e is not to be confounded with the coupling 
constant e defined earlier). Eq. (P^D indicates that at a certain value of the mean 
density of obstacles ^Jq the mean field solution will break down. However, the specific 
value obtained for the upper bound is a result of the mean field expression for the 
critical area Ac, Eq. (pT]) , and should not be taken as a quantitative, but a qualitative 
result. Nechaev et al. have stressed that the mean field approximation is valid in 
the vicinity of the boundary curve defined by the lower bound inequality Eq.(p9D 
1 5^ (see FIG.(l)). The present approach using the area of the loop as the order 
parameter and giving the new upper bound gives strong support of this result that 
the mean field solution is restricted to a small neigborhood above the boundary 
curve. 
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IV Final result and discussion 



It has been argued in that the collapsed phase can be identified with a randomly 
branched polymer. Here we give explicit support of this idea. Consider Eq. (|27| ) and 



substitute p = ^ in the mean field approximation where > N^., i.e. above the 
critical length. We obtain: 

A = ^N(l-^ln{N) + ..}j. (33) 
Taking the terms in parentheses as the first powers of an exponential, one finds: 

A = —N{N~^2) = (34) 

Exploiting the result of Cardy ||29| for SAWs in = 2 that {A) ~ (-R^), one obtains 
for i? = ((i?2))l: 

R ~ IN^ (35) 

This is the scaling behavior for randomly branched polymers without excluded vol- 
ume interaction. It is therefore very likely that the collapsed phase of the loop 
corresponds to a randomly branched polymer. 

Note that the result of Eq.(0) corresponds to the free part of the free energy 
in Eq.(^). i.e. without excluded volume. This is valid as a first approximation 
because the area A is not directly coupled to the density p = |0p at the level of 
the free energy, and the excluded volume is not renormalized by the topological 
interactions in the mean field approximation. Moreover, a Landau expansion in 
terms of segment density p and tangent vector density variables indicates that 
the p and decouple at first order because /c^j^(k) = 0, while interactions occur 



only at higher order. |14| 

As a consequence, the two-dimensional excluded volume must still be taken into 
account. That yields the well known d = 2 branched polymer scaling [^6|,^ for the 
area or the mean sqare end-to-end vector respectively, i.e. ~ IN^. 



V Outlook to the three-dimensional problem 

Finally, let us point out some possible future perspectives for the d = 3 entanglement 
problem by coming to the mathematical difficulty of finding a correct knot invariant. 
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While the known knot invariants in their algebraic form seem to be only of limited 
use in the polymer context, the work of Witten [|ll|] showing an equivalence of the 
Jones polynomial (actually a more general one) and (in general non-abelian) Chern- 
Simons field theory has brought the knot problem closer to physics again. Witten 
showed that the expectation value of Wilson lines averaged over a Chern-Simons 
action functional integral gives a knot invariant for framed links. This invariant can 
also be reproduced perturbatively giving the sqare of the Gauss invariant as a first 
approximation and higher order knot invariants |^||31]] ||12|| . 



In addition it has been shown recently that the invariants appearing in Chern- 
Simons perturbation theory are intimately related to so-called Vassiliev invariants 
(see e.g. ||3^ and the literature quoted therein). 



These results may open new ways of solving the polymer entanglement problem 
by a "topological perturbation theory". It might give a range of validity for using 
the Gauss invariant for ensemble of random walk chains or rings. 
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FIG. 1: The phase diagram of the collapse transition. The shaded area bounded 
by the sohd curve corresponds to the collapsed state of the loop as obtained in [22] . 
The new upper bound in Eq.(32) gives rise to the new boundary curve (dashed line). 
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